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Abstract 

We study long-range percolation on the hierarchical lattice of order N, where any 
edge of length k is present with probability = 1 — exp(— f3~ k a), independently of 
all other edges. For fixed /3, we show that the critical value a c ((3) is non-trivial if and 
only if N < (3 < N 2 . Furthermore, we show uniqueness of the infinite component and 
continuity of the percolation probability and of a c (/3) as a function of f3. This means 
that the phase diagram of this model is well understood. 

1 Introduction and main results 

The use of percolation theory in statistical physics has long been recognized. The study of 
long-range percolation on 7L d goes back to [18] and led to a series of interesting problems 
and results [21 [T7J HJ El [19] ; see [6j Section 2] for an extensive overview. In [10J asymptotic 
long-range percolation is studied on the hierarchical lattice (to be denned below) for 
iV — > oo. The contact process on fl^ for fixed N has been studied in [3]. 

In this paper we study the case of finite N. Long range percolation on the hierarchical 
lattice is quite different from the usual lattice: classical methods break down and results are 
different. The results and methods in this paper could appeal to both mathematicians and 
physicists. 

For an integer N > 2, we define the set 

tt N := {x = (xi, x 2 , . . .) : Xi e {0, 1, . . . , N - 1}, x t < oo}, 
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Figure 1: Hierarchical lattice (the "leaves") of order 2 with the metric generating tree at- 
tached. 



and define a metric on it by 

'O ifx = y, 

max{i : x { ^ y { } if x ^ y. 

The pair (Q N ,d) is called the hierarchical lattice of order N. 



d(x,y) 



One can think of the vertices in the hierarchical lattice as the leaves of a regular tree without 
a root, see Figure HJ The metric d can then be interpreted as the number of generations 
(levels) till the "most recent common ancestor" of two vertices. Let N be the non-negative 
integers, including and M + := N \ {0}. The set Qn is countable, and we can introduce a 
natural labeling of its vertices via the map / : ~ ► N given by 

oo 

i=i 

We will sometimes abuse notation and write n for / _1 (n) G Qn- 

The metric space (fljv,d) satisfies the strengthened version of the triangle inequality 

d(x,y) < max(d(x, z), d(z, y)) 

for any triple x, y, z e Q^. Such spaces are called ultrametric (sometimes non- Archimedean) 
[20] . For i 6 fljv, define B r (x) to be the ball of radius r around x. Several important 
geometrical properties follow from the definition of the space (£In, d) and its ultrametricity: 
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1. B r (x) contains N r vertices for any x; 

2. for every x G Qn there are (N — l)A^ fe_1 vertices at distance k from it; 

3. if y G B r (x) then B r (x) = B r (y); 

4. as a consequence of the previous property, for all x, y and r we either have B r (x) = 
B r {y) or B r (x) n B r {y) = 0. 

Now consider long-range percolation on Q^. Every pair of vertices (x,y) G Qn x &n is 
(independently of all other edges) connected by a single edge with probability 



where k = d(x, y) and where < a < oo and < (3 < oo are the parameters of the model. 
The vertices x G and y G Qn are in the same connected component if there exists a 
path from x to y, that is, if there exists a finite sequence x — Xq, x±, . . . , x n = y of vertices 
such that every pair of points with 1 < i < n shares an edge. The edges are not 

directed. 

We denote the size of a set S of vertices by \S\. The connected component (also called 
"cluster") containing the vertex x is denoted by C(x). Since |C(x)| has the same distribution 
for every x G Qn we may study |C(0)| instead of |C(x)|. 

Let Pq,^ be the probability measure governing this percolation process (on the appropri- 
ate probability space and sigma- algebra) and the corresponding expectation operator. 
When no confusion is possible, we omit the subscripts a and j3. Denote 



It follows from a standard coupling argument that 6(a, (5) is non-decreasing in a for any 
given (3. Therefore, it is reasonable to define 



Throughout the paper we use the following notation. For a set S of vertices, let S := 
\ S denote its complement. The set C n (x) is the cluster of vertices that are connected 
to the origin by a path that uses only vertices inside B n (x). For disjoint sets Si, S 2 C Qn, 
the event that at least one edge connects a vertex in Si with a vertex in S2 is denoted by 
Si <H> S2. The notation Si ft S2 denotes the event that such an edge does not exist. Let 




0(a,/9):=P ai/J (|C(O)| = oo). 



(1) 



a c (/3) := inf{a > : 9(a,/3) > 0}. 
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C™(x) be the largest cluster in B n (x); if more such clusters exist, C™(x) is defined to be one 
of them, chosen uniformly among all possible candidates. In any case, 

IC(*)I= \C n {y)\. 

Theorem 1 ((non) -triviality of the phase transition) 

(a) a c (P) = 0for/3<N, 

(b) < a c (/3) < oo for N < f3 < N 2 , 

(c) a c ((3) = oo for > N 2 . 

Theorem 2 (uniqueness of the infinite cluster) There is a.s. at most one infinite cluster, 
for any value of a and j3. 

Theorem 3 (continuity of 6) The percolation function 9(a,/3) is continuous whenever a > 
0. 

Theorem 4 (continuity ofa c ((3)) The critical value a c ((3) is continuous for G (0,iV 2 ) and 
strictly increasing for j3 G [N, N 2 ). Finally, a c ((3) /* oo for j3 /* N 2 . 

In order to prove Theorems [3] and HJ we need the following result, which is interesting in 
its own right. 

Theorem 5 (size of large components) If a and (3 are such that 9 := 9(a,(3) > 0, then for 
every e > 0, 

hmP Q ,,(|C(O)|>(0-^) = l. ( 2 ) 

In the next two sections we prove Theorem [1] and Theorem [2] respectively. After that we 
prove the remaining results, and we end with a discussion about possible generalisations. 



2 Proof of Theorem H 

Proof of (a). Denote by Ek the event that the origin shares an edge with at least one vertex 
at distance k. Then 

P(£ fc ) = l-exp f-^(N-l)N k - 

and the events (Ek)k>i are independent. It is easy to see that if /3 < N then Ylh=i ^(^k) 
diverges for any a > 0. Therefore, by the second Borel-Cantelli lemma, infinitely many of 



4 



the events E k occur with probability 1 and 9 (a, (3) = 1, for any a > and < (3 < N. This 
implies a c {(3) = for < (3 < N. □ 

Proof of (c). By monotonicity it suffices to prove that a c (N 2 ) = oo, so we now take (3 = N 2 . 
A straightforward computation shows that for every j, 



i+k-l 



m(v) » Bj(v)) = 1 - exp ^-a^L^i £ ^^ryj = 1 - exp (-£) , 

which is strictly less than 1, and independent of j. Let n = and n i+1 = inf{n > rii : 
B ni (0)^^(0)}. Since 



{C(0) = oo}cn- {^(0)oB ni (0)}, 



it is enough to prove that there a.s. exists i such that £> n ^(0) f4 -B ni (0). Because the events 
{B ni (0) -H- <B„ t (0)} are independent and all have the same probability strictly less than 1, 
the result follows. □ 

Proof of (b). The strict positivity of a c ((3) follows from the fact that 

OO OO , \ 

^(iV-l)iV fc -V. = E(^-l)^- 1 l-exp(-^) < 
k=i k=i ^ p ' 

^ a(N-l)^ fN\ k /Ar 1S 1 

which can be made strictly smaller than 1 by choosing a small enough. Hence the expected 
number of edges from a given vertex is strictly smaller than 1, and by coupling with a 
subcritical branching process, the almost sure finiteness of the percolation cluster follows. 

The second inequality is much more involved. Choose an integer K and a real number r\ 
such that 

y^< V <(N K -l) 1/K ; (3) 

this is possible since \f(3 < N. We say that a ball of radius nK is good if its largest connected 
component has size at least i] nK . Denote by s n the probability that a ball of radius nK is 
good, that is, 

Sn :=P(&(0)|>^)- 

By convention, we set s = 1. 

We say that a ball of radius nK is very good if it is good and in addition its largest 
component shares an edge with the largest component of the first (the one with the smallest 
index) good sub-ball in the same ball of radius (n + 1)K. Note that according to this 
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definition, the first good sub-ball of diameter nK in a ball of radius (n+l)K is automatically 
very good. 

Since (N — 1) > rj K , Bf n +i)K:{Q) w ih certainly be good if 

(a) it contains N K — 1 good sub-balls of radius nK, and 

(b) all these good sub-balls are very good. 

We next estimate the probability of the events in (a) and (b). 

Clearly, the number of good sub-balls of radius nK in a ball of radius (n + 1)K has a 
binomial distribution with parameters N K and s n . Furthermore, given the collection of good 
sub-balls, the probability that the first such good sub-ball is very good is equal to 1, and the 
probability for any of the other good sub-balls to be very good is at least 



since the distance between two vertices in a ball of radius (n + 1)K is at most (n + 1)K and 
the largest component of a good sub-ball contains at least n nK vertices. We conclude that 
the number of very good sub-balls is stochastically larger than a random variable having a 
binomial distribution with parameters N K and s n (l — e n ). It follows that 



where Bin(n,p) denotes a random variable with a binomial distribution with parameters n 
and p. Notice that since 1/x > e~ x for all x > 0, we have 




and hence for any S > we can find a so large that e n < 5 n for all n. 

The expression in (j3J) is very close to the iteration formulae in [11], and it pays to recall 
the setup from that paper first. 

Let, for some given m, the map 





(4) 



7r(p) = F(Bin(m,p) > m — 1) 



be defined. Define Uq — 1 and, for neN 



u n+ i = -n{pu n ). 
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Writing G p (-) for n(p-) we then obtain 

u n+1 = Gp{u n ). (5) 

In [II] . this iteration arises in the study of fractal percolation and it is shown that the limit 
u = lim^oo u n always exists and is positive if and only if p is so large that the equation 
G p (x) = x has a positive solution. This result is very similar (and can be proved in the same 
way) as the classical non-extinction criterium for branching procesess. When p approaches 
1, then the largest solution of G p (x) = x also approaches 1, at essentially the same rate, and 
therefore also the limit u approaches 1. 

Now observe that (j3J) can be rewritten as 

Sn+l > Gt-e^Sn). (6) 

This is very similar to (151) , the only difference being that the subscript of the iteration 
function depends on n now. However, we already showed above that for a large enough, e n 
goes down exponentially fast at any given rate. It is not hard to believe that this implies 
that s n converges to 1 exponentially fast, and we make this precise now. 

We have 

F(Bin(n,p) > n - 1) > 1 - 
and writing C = ( N 2 ) we arrive at the inequality 

s n +i > 1 - C(l - s n + s n e n f > 1 

Writing £ n = 1 — s n this gives 

in+i<C^ n + e n f. (8) 

Choose first 7 so small that 4C < 7~ x and then a so large that e n < 7™ and £1 < j 2 . In an 
inductive fashion, if < ^ n+l then 

£n+i < C(£ n + e n f < 4C( 7 n+1 ) 2 < 7 2n+1 < 7 n+2 , (9) 

which implies that £ n < 7™ +1 for all n. Hence, for a large enough, s n converges to 1 
exponentially fast. 

The exponential convergence of s n to 1 is not quite enough for our purposes. Indeed, s n 
represents the probability that a ball of radius nK contains a component of size at least r\ nK , 
but this component does not necessarily contain the origin. Therefore, we have to make one 
extra step. Let 

t n :=F(\C nK (0)\> V nK ). 



n 



2 ,a-pr, 



C(l - s n + e r 



(7) 
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We claim that 

tn+i >t n x ¥{Bm{N K - 1, s n (l - e n )) > N K - 2). (10) 

To see this, we argue as before. If |C n ^-(0)| > r] nK , then B n K(0) will be the first good sub- 
ball in the derivation above. If this component is connected to at least N K — 2 other large 
components in £>( n+ i)£-(0) as above, then the component of the origin in £>( n+1 )^ is large 
enough, that is, has size at least ^( n+1 ) A ". From this, ( TTUj) follows. Since a simple coupling 
gives that 

F(Bin(N K - 1, s n (l - e n )) > N K - 2) > F(Bin(N K , s„(l - e n )) > N K - 1), 

and since the derivation above actually gives that the right hand side of this inequality 
converges to 1 exponentially fast, it follows that 

lim t n > 0, 

n— ¥oo 

which is enough to prove the result. □ 

Remarks (I) By the proof of the strict positivity of a c (/3) for > N, we also may deduce 
that a c (P) is not differentiable in N. Indeed, a c (/3) = for < N and since 

oo 

J2(N - l)iV fc -V < a(N - 1)— — 

k=l P 

for (3 > N, we have a c ((3) > for all (3 > N. This in turn implies that for all /3 > N we 
have 

q c Q9) - a c (N) 1 
p- N ~ N -I 

(II) Since we may choose 7 arbitrary small in equation we in fact have that for every 

e > 0, we can choose a so large, that t n > 1 — e. This implies that for every G (N, N 2 ), 
we can choose a so large such that 9(a, (3) > 1 — e. 



3 Proof of Theorem [2 

We will use Theorem from [IB] : 

Theorem 6 (Gandolfi et al. Jig)/) Consider long range percolation on 7L d with the properties 

1. the model is translation-invariant, and 

2. the model satisfies the positive finite energy condition. 
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Then there can be a.s. at most one infinite component. 



In order to be able to use this result, we will first embed the metric generating tree into 
Z in a stationary (and ergodic) way. The embedding will be such that for each r, we have 

a. any ball of radius r will be represented by N r consecutive integers, 

b. the collection of balls of radius r partitions Z. 

We first describe the construction rather loosely, and after that provide a formal construction. 
For ease of description, a collection of m consecutive integers is called an interval of length 
m. 

The ball of radius 1 containing 0, that is, £>i(0) is chosen uniformly at random among all 
N possible intervals of length N containing the origin of Z. Once we have chosen this ball, 
all other balls of radius 1 are determined by requirements (a) and (b) above, although it is 
not yet clear at this point to exactly which balls in they correspond. To get an idea of 
this first step of the procedure, note that for N = 2 there are only two possibilities, one of 
which is depicted here: 



The other possibility is obtained by translating the edges over one unit to the right (or to 
the left, for that matter). 

Next, we determine £> 2 (0). The ball £> 2 (0) is a union of N balls of radius 1 and contains 
Bi(0). There are N possible ways to achieve this, keeping in mind that any ball of radius 2 
must - according to (a) above - be an interval of length N 2 . We now simply choose one of 
the N possible ways to do this, with probability 1/N each. Once we have chosen £> 2 (0), all 
other balls of radius 2 are determined for the same reason as before. The following picture 
illustrates a possible choice for B 2 (0) given the choice of B±(0) made before. 



We can continue this procedure as long as we wish, and in doing so we obtain a metric 
generating tree which is isomorphic to the tree depicted in Figure 1. This last statement 



s\ s\ s\ 



-3-2-101234 




-3-2-101234 
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perhaps requires some reflection: one can see that this holds by first identifying the two O's 
in both graphs, and then build up the balls B r (0), r — 1, 2, . . ., in that order. 

It is intuitively clear that this construction yields a stationary metric generating tree in 
the sense that the distribution of the stochastic process which assigns to each pair {z, z'} of 
points in Z the distance between them, is invariant under integer translations. However, we 
would like to formalise the construction in such a way that not only stationarity follows as 
an easy corollary, but we also obtain that the embedding of the metric generating tree is in 
fact ergodic with respect to translations. 

A possible formal construction is the following. Our probability space is the unit interval 
[0,1] endowed with Lebesgue measure on its Borel sigma field. For 7 e [0,1], let 7 = 
O.7172 ■ ■ • be its iV-adic expansion, that is, 



where we ignore those 7 for which the expansion is not unique - this is a set of Lebesgue 
measure zero anyway. In the construction above, we saw that for each r, S r _i(0) can be seen 
as one of the balls of radius r — 1 among the balls making up B r (0). The metric generating 
tree corresponding to 7 G [0, 1] is obtained as follows. We let £> r (0) be such that <B r _i(0) is 
the (j r + l)-st ball in £> r (0), counted from left to right. For instance, in the preceding two 
figures with N = 2, we have that 71 = 1 and 72 = 0. The map which assigns to each (apart 
from the exceptional null set discussed before) 7 a metric generating tree is denoted by 0. 
This map is invertible on a set of full Lebesgue measure. 

It is clear that this construction formalises the informal description given earlier. Further- 
more, one can write down explicitly the transformation S : [0, 1] — > [0, 1] which corresponds 
with the left-shift T on the space of metric generating trees in the sense that <p o S = T o 0, 
hence T = <pS(j)~ l . Indeed, a little reflection shows that S can be described as follows: if 
Y(j) = min{n; 7„ 7^ iV — 1} then (that is, the k-th digit in S(j)) is given by 



This transformation has been studied in the literature and goes by the name Kakutani - 
Von Neumann transformation, see e.g. [12]. It is easy to check that Lebesgue measure is 
invariant under the action of S, and this immediately proves that the construction of our 
random metric generating tree is stationary on Z. 



oc 



7 = ^7niV-, 



71=1 




iffc<F( 7 ), 
7fc + l iffc = F( 7 ), 
7 fc if k> Y{i). 
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Proof of Theorem^ The construction above shows that the metric generating tree can be 
embedded into Z in a stationary way. We claim that this implies that the whole long range 
percolation process on the hierarchical lattice can be realised as a stationary percolation 
process on Z. To see this, we assign a uniformly- [0, 1] distributed random variable U e to 
each edge e in such a way that the collection is independent. 

Given a realisation of the metric generating tree, we declare edge e to be open if U e < 
1 — exp(— ct//3' e '), where |e| denotes the length of e. This gives a realisation of the percolation 
process with the correct distribution, and shows that we have embedded the full percolation 
process on the hierarchical tree in a stationary way. Since every pair of vertices shares an 
edge, with positive probability, irrespective of the presence or absence of other edges, the 
positive finite energy condition is met and the result follows. □ 

With a little more work one can also see that the construction is in fact ergodic, that is, any 
event which is invariant under the shift on Z has probability or 1. To show this, we first 
show that Lebesgue measure on [0, 1] is ergodic with respect to S. This result is known, but 
we give a simple (and new) proof for the convenience of the reader. 

Lemma 7 Lebesgue measure on [0, 1] is ergodic with respect to S . 

Proof. Consider the first digit in each of 7, £(7), S 2 ^), • • •• From the construction we imme- 
diately conclude that the first digit follows the periodic pattern 0, 1, ... , N—l, 0, 1, 2, ... , N — 
1, . . . starting at any number. Hence the first digit is just adding 1 modulo N. The second 
digit can only change when the first digit is an N — 1 , and then it also changes according to 
adding 1 modulo N. In general the k-th digit can only change when the (k — l)-st digit is an 
N—l, and then the change consists of adding 1 modulo N. It follows from these observations 
that the orbit of 7 under the action of S visits any iV-adic interval I m ^ = [kN~ m , (k+l)/N~ m ] 
with frequency N~ m , that is, 



where 1a denotes the indicator function of A. 

Now consider the collection Ai of invariant probability measures for the transformation 
S. From the fact that Lebesgue measure preserves measure under S we see that Ai is not 
empty. It is well known and easy to see that the set M. is convex, and that the ergodic 
measures are precisely the extremal points of M.. Since Ai 7^ 0, this implies that there is at 
least one ergodic measure with respect to S. 




n-l 



(11) 



3=0 
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Let v be any ergodic measure with respect to S. It then follows from the ergodic theorem 
and (JTT1) that v[l m ^) = N~ m for any m and k — 0, . . . , N m — 1. However, there is only one 
measure that satisfies this condition, namely Lebesgue measure on [0, 1]. Hence v must be 
Lebesgue measure, which we already know is indeed invariant. □ 

Theorem 8 The embedding of our long range percolation process on the hierarchical lattice 
into Z is ergodic. 



Proof. From Lemma [7] it follows that the metric generating tree is embedded ergodically. 
Adding the i.i.d. random variables U e as before does not destroy ergodicity, and the final 
configuration is a factor of this ergodic process and hence ergodic itself. □ 



4 Proof of Theorem [5 

The proof consists of three steps: 

1. For every constant K > the indicator function of the event that both |C(0)| = oo 
and |C n (0)| < K(f3/N) n converges a.s. to as n — > oo. 

2. The fraction of the vertices in B n (0) which are in a cluster of size at least K((3/N) n , 
converges a.s. to 9 as n — > oo. 

3. Combine the previous two steps. 



Step 1. We compute 



P ( |C(0)| = oo \{n e N; |C n (0)| <K[^) }| = oo 



Let ni be the smallest n for which |C„(0)| < K(f3/N) n , if C ni (0) B ni (0), then n i+1 is the 
smallest n>m for which C^(0) B n (0) and for which |C n (0)| < K((3/N) n . 
Since |C ni (0)| < K((3/N) n \ we have 

p(c ni (0)^B^(0)) < p(Ci(0)^^(0)||Ci(0)| = [k(§) 

( 00 T 

< l-e XV (-aK(§yj2(N-l)—) (12) 



Ni- 1 



1 — exp (^—aK 



j=i+i 



N-l 
0-N 
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The right hand side is strictly less than 1 and is independent of rij. So there will be an 
for which {C„ t (0) ft B ni (0)}, and it follows that 

P (jC(0)| = oo\\{n e N; |C n (0)| < K (JA }| = oo) = 0. 

Step 2. We use the random embedding of the hierarchical lattice in Z, introduced in the 
previous section. By Theorem [H] and the ergodic theorem we have, for every k > 0, 



2N n + 



— y 

- + 1 ^ 




a.s. 

y 



r (r\[w°)\> K (ji) 



\j=k 

From step 1 we know that this final probability increases to 9 as k — > oo, and it follows that 
A(n) := Y l( \Cn(x)\ > K f 41 1 ^• 

x=-N n v \ / / 

Note that the collection vertices {— N n , —N n + 1, —N n + 2 • • • , N n } contains the image under 
the embedding of the ball B n (0) and this image contains a fraction N n /(2N n + 1) of those 
vertices. Whether or not |C n (a;)| > K(fi/N) n is independent for vertices in different n-balls, 



so 



and A 2 (n) := — A\{n) are independent. Furthermore, Ax(n) and A 2 (n) are bounded 
above by 1 and have asymptotically the same mean. Combining these observations with 

A l (n)+A 2 (n) = A{n) ^-^6 

n—too 

finishes the proof of Step 2. 

Step 3. The strategy is to split those components in £> n+1 (0) which are at least of size 
K(f3/N) n into clusters roughly of size K((3/N) n . Then we use those clusters as "meta- 
vertices" for an iV-partite graph, in which meta-vertices in different n-balls are connected if 
the clusters they represent are connected by an edge of length n + 1. Meta-vertices in the 
same n-ball never share an edge. We show that if we choose K and n large enough, then the 
largest component of the graph of meta-vertices contains a fraction of the meta-vertices close 
to 1, which shows that for large n, the fraction of vertices in the largest cluster of i3 n+ i(0) is 
close to 9. We will be more precise now. 
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By step 2 we know that for every K > 0, every e > and all large enough n, it holds 
that 



We condition on the event that all n-balls in B n+ i (0) are good. The probability of this event 
is bounded below by (1 — e) N > 1 — Ne. 

Now, for every good ball B n (y), y G Qn, we make a partition of the set 



in "met a- vert ices" . For the moment we denote this set by B' n (y). For x G B' n (y) we make 
a partition of C n (x) in [\C n (x)\ / (\ K (/3 / N) n ]) \ sets, which all have size at least \K((3/N) n ] . 
Here \x] := inf{n G Z; n > x} is the ceiling of x and |_^J : = sup{n 6 Z; n < i} is the 
floor of x. The vertices that are not in such a cluster are ignored for the moment. Denote 
the collection of meta- vertices that contain vertices in B n+1 (0) by V n . We note that if 
B n (y) is good and K is large enough, then it contains at least (6 — e)N n /\2K(/3/N) n ] > 
(6 - e)N n /(3K(/3/N) n ) vertices. 

We construct a new iV-partite graph on V n as follows. Let V n be the vertex set and let 
£ n be the set of edges between those vertices. This edge set is obtained as follows. Choose 
\K(/3/N) n ] original vertices from every meta-vertex V n . Chosing those vertices may be done 
in any way that is independent of the presence of edges of lenght n + 1 or larger. Denote 
these sets by A n - The meta-vertices x, y G V„ share an edge in £„, if there is at least 1 edge 
in the original graph that is shared by vertices that make up the sets in A n corresponding 
to x and y, and if the original vertices that make up x and y are at distance n + 1 of each 
other. Otherwise there is no edge between the meta-vertices. 

As observed before, the number of meta-vertices in V„ that consist of vertices from a 
good ball B n (x), is at least (9 - e)N n /(3K(/3/N) n ). Since f3 < N 2 , this quantity grows to 
oo as n — t> oo. The expected degree of a vertex in V n exceeds 




We now fix e. The ball B n (y) is said to be good if 






which is larger than A := (N — l)(0 — e)aK/(6(3), for all large enough n. This holds for every 
K > 0, and therefore the expected degree can be chosen to be arbitrary large. 
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This iV-partite graph falls within the class of inhomogeneous random graphs of Bollobas, 
Janson and Riordan [7]. The degree of every meta- vertex is asymptotically Poisson dis- 
tributed, with mean bounded below by A and we know that the (unique) largest component 
of such an iV-partite graph contains with high probability (in the limit for n — > oo) a fraction 
p of the met a- vert ices, where p is the largest solution of 

1 - p = e~ xp . 

By tuning K, A can be chosen arbitrary large and p can be taken such that p > 1 — e. So, 
for every e > and large enough n the graph (V n ,£ n ) contains a unique giant component, 
which contains a fraction (1 — e)N of the vertices in V n , with probability at least 1 — e. 

Since we have conditioned on the event that all n-balls in £> n+1 (0) are good, the fraction 
of vertices in £> n+ i(0), that are part of vertices in V n is bounded below by 9 — 2e. (The 
factor 2, is due to the fact that the sizes of different meta- vertices differ at most a factor 2). 
Therefore, with the same conditioning, the largest cluster in $ n+1 (0) is at least of size 

(p - e){9 - 2e)N n > (1 - 2e)(0 - 2e)N n 

with probability exceeding 1 — e. Now multiplying by the probability that all n-balls in 
B n +i(0) are good, gives that the probability that the largest cluster in B n +i(0) is at least of 
size (l-2e)(9-2e)N n is bounded below by (l-e)(l- JVe). By chosing e' G (0, ej max(4, N+ 
1)), we obtain that P(|C n (0)| > (9 — e')N n ) is at least 1 — e' and this finishes the proof. □ 

Remark We realize that it is possible to prove the statement of Step 2 by using the strong 
law of large numbers. If we do this, then it is only a small step from the proof of Theorem [5] 
to a proof of Theorem [2j However, we think that the proof presented in the previous section 
contains some valuable ideas and therefore should be included in this paper. 



5 Proof of Theorem [3 

Continuity proofs of percolation functions typically split into separate proofs for left and right 
continuity, one of which typically follows from standard arguments. In this case, continuity 
from the right in a and continuity from the left in /3 are the easy parts: 

Lemma 9 9 (a, (3) is continuous from the right in a > and continuous from the left in 
/3 > 0. 

Proof. For a > and (3 < N, 9(a,(3) = 1, so the statement of the lemma holds in that 
domain. Note that 



{d(0) <-> Bi(0)} C {^(0) ^ Bi_i(0)} 
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and 



{|C(0)| 



00} = nZ {Ci(o) Bi(o)}- 



(13) 



A straightforward computation yields that for f3 > N, 




N- 1 fN 




) 



(14) 



Since |Cj(0)| depends on the state of only finitely many edges, this expectation is continuous 
in a and (3. In particular, P(Cj(0) <-> &i(0)) is continuous from the left in (3 for (3 > N 
and therefore it is continuous from the left for /3 > (in fact, it is also continuous from the 
right). Furthermore, P(Cj(0) -H- £>«(0)) is increasing in a, decreasing in (3 and decreasing in i. 
Because a decreasing limit of increasing (resp. decreasing) functions, which are continuous 
from the right (resp. left) is continuous from the right (resp. left), the statement follows. □ 

In order to prove that 6(a,f3) is continuous from the left in a > and continuous from 
the right in (3 > 0, we use a renormalisation argument. Fix a > and N < (3 < N 2 . To get 
insight in the argument we first (falsely) assume that for given e > and large enough finite 
k, there is a 5 > such that, 



(Although the assumption is false, we can get this probability arbitrary close to 1, by choosing 
k large enough, 6 small enough and using Theorem [51 ) 

Now we use renormalisation. The balls of radius k are considered as vertices of which 
we call "met a- vert ices" . If two vertices in the original model have distance k + I, then the 
meta-vertices in which they are contained are at distance I. Vertices in the new model are 
connected if and only if the largest clusters in the original fc-balls, represented by these 
vertices, are connected by an edge. The new model is again a percolation model on fl^. 

Let x and y be meta-vertices, at distance / of each other. Define, for 5 > small, 



Given the states of all other edges, the (conditional) probability that x and y are connected 
to each other is always bounded below by 



P Q _ w (|Cno)| > (0(a,P) -e)N k ) = 1. 



(a-5)((6(a,{3)-e)) 2 



N 2k 



a := 



(f3 + 8) k ' 



1 - exp(-(a - 5)((8(a, (3) - e)N k ) 2 {(3 + 5)~ {k+l) ) 



and by the choice of a', this is just 



1 -exp(-a'(/3 + 5)-'). 
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Hence, the renormalized model stochastically dominates the percolation model with param- 
eters a' and (3 + 5. 

Since N 2 / (0 + 5) > 1, a' can be chosen arbitrary large by choosing k large. In particular 
it can be chosen such that 8(a',f3 + 5) > 1 — e, (by the second remark after the proof of 
Theorem [1]) . It follows that for large enough k 

V«s,p+s (\C\ = oo) > 6(a', (3 + <T)P a _ w (0 £ C fc m (0)) 

> (l-e)($(a,l3)-e) 

> 9(a,/3)-2e. 

The only problem is that we have incorrectly assumed that 

V a -8,fi + 6 {K(0)\ > (6(u,(3)-e)N k ) = 1, 

and we will deal with this problem now. We need the notion of mixed percolation (cf. [8]). 
Mixed percolation involves independently removing vertices, together with all of its adjacent 
edges. Formally, the measure P™^^ is constructed as follows. A vertex in is open with 
probability 1 — 7, independently of the states (open or closed) of the other vertices in Q^- If 
x and y are both open vertices then they share an edge with probability 1 — exp(— a//3 d ( x ' y >), 
independently of which other edges are present. 

Lemma 10 Let /3 > N. For every e > 0, there exists 7 > ; such that 

p^(|C(o)| = 00) < ps8S )A7 (|C(o)| = 00). 

Before giving the proof of this result, we show how it can be used to prove Theorem [3j 
The following lemma suffices. 

Lemma 11 If 6(a,(3) > 0, then for every e £ (0,9(a, j3)), there exists 5 > such that 
9(a-5,/3 + 5) > 9(a,fi)-e. 

Proof. Fix a, (3 and e > 0. Let a' be such that 

9(a', (N 2 + f3)/2) > 1 - e/3, 

which is possible by Theorem [H(b) and the second remark after its proof. Furthermore, let 
7 £ (0, e/3) be such that 

9(a', (N 2 + f3)/2) < P^ 2+/3)/2i7 (|C(0)| = 00), 

which is possible by Lemma [10J Let K be such that the following conditions are satisfied: 
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1. a(6(a,P)-e/2)2(N 2 /P) K >3a', 

2. P o ^(|Ck(0)| > (9(a, 0) - e/3)N K ) > 1 - T /2, 

which are possible by respectively iV 2 > and Theorem [51 Finally, let 5 > be such that 
5 < min(a/3, (N 2 - 0)/2) and 

Pa-W|C£(0)| > (9(a,0)-e/3)N K ) > 1 - 7j (15) 

which is possible by the continuity of the probability in a and for finite if. 

We say that the ball Bk{x) is good if has size at least (9(a,f3) — e/3)N K . Delete 

all vertices that are in a ball of diameter K which is not good and also all vertices that are 
not in the largest cluster of good balls. As above, we interpret the remaining components as 
the vertices of the hierarchical lattice of order A" in which vertices are independently deleted 
with probability at most 7, by (fl5|) . Remaining clusters in the original graph, of which the 
vertices are at distance K + I, are connected by at least one edge with probability at least 

1 - exp(-(a - S)(9(a, 0) - e/3)2N 2K P~^ K+l) ) > 1 - exp(-2a'p~ l ), 

irrespective of the existence or absence of other connections. Here we have used that a — 5 > 
2a /3. Hence the rescaled process stochastically dominates a mixed percolation process with 
parameters 2a', (3 and 7. 

Now note that by exchangability 

F a _ Si p +s (\C K (Q)\ > (6(a,f3)-i/3)N K ) 

> (e(a,P)-e/3)F a - S ,f, +s (\C%{0)\ > (0(a,P) - e/3)N K ) 

> (0(a,/3)-e73)(l- 7 ). 

Furthermore, conditioned on being in the largest cluster of a good ball, the probability 
that is in an infinite cluster if the parameters are a — S and ft + 5 is larger than 1 — e/3. 
Combining these observations and 7 < e/3 gives 

9(a -5,0 + 8)> (9(a, 0) - e/3)(l - e/3)2 > 9(a, 0) - e. 

□ 

It remains to prove Lemma (TOj Before giving a proof of this lemma we define a directed long- 
range percolation model and relate this to the undirected model. In the directed version, 
vertices in Qn are open with probability 1 — 7. If vertex x is open, then a directed edge from 
x to y is present with probability 1 — exp (— a(3~ d<yX ' y ^ . Conditioned on the states of the 
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vertices (open or closed) the presence or absence of an edge is independent of the presence or 
absence of other edges. The corresponding measure we denote by P™^ d . The set of vertices 
which can be reached by a path from vertex x is denoted by C(x). Note that in th directed 
model, the presence of a path from x to y does not necessarily imply that there exists a path 
from y to x. We define the directed version of the original (not mixed) measure, P a)/ g, in a 
similar way and note that F a ^ = P™^^- Standard arguments (see e.g. [HI EES]) can be used 
to show that 

K^ d (\C(0)\ = oo) = Kj° d (\C(0)\ = oo). 

Proof of Lemma\TQ The directed mixed percolation graph with parameters a, /3 and 7 can 
be obtained as follows (the ordinary model can be obtained by taking 7 = 0). We assign 
i.i.d. random variables X x to the vertices 1 6 Hjv, all Poisson distributed with parameter 
a(N — l)/(/3 — N). We construct a directed multi-graph (a graph in which multiple edges 
between two vertices in the same direction are allowed). Vertices are open with probability 
1 — 7, independently of each other. If x is open, then X x directed edges start at x. The 
endpoints of these edges are independently chosen from Qn \x, and a vertex at distance r 
of x is chosen with probability (/3 — N)(N — l)~ l [3~ r . If x is closed, then no edges start at 
x. We obtain the original directed graph by replacing the collection of all edges from x to y 
(if there is at least one) by a single edge from x to y, for all x,y £ Oat- 

Let Z\ be a Poisson distributed random variable with parameter a(N — l)/((3 — N). 
Furthermore, let Z 2 = Y1Y2, where Y\ is equal to 1 with probability 1 — 7 and equal to with 
probability 7, and where Y 2 is independent of Y x and Poisson distributed with parameter 
a(l + e)(iV — 1)/(/3 — N). For the ordinary percolation model the number of edges starting at 
x in the multigraph is distributed as Zi, while for the mixed percolation model, the number 
of edges starting at x is distributed as Z 2 . It is now easy to check that for e > there is a 
7 > 0, such that ¥(Zi = 0) = P(Z 2 = 0) and for this 7 and all k > we have, 

F(Z 2 > k\Z 2 > 0) = F(Y 2 > k\Y 2 > 0) > P(Zi > k\Z x > 0). 

The statement of Lemma [TU1 now follows by a straightforward coupling argument. □ 



6 Proof of Theorem [4 

The proof of Theorem H] is split into separate proofs of continuity from the right and from 
the left of a c (/3). 

Lemma 12 a c ((3) is strictly increasing on (3 £ (N, N 2 ) and continuous from the right on 
/3£(0,iV 2 ). 
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Proof Theorem [3] implies that 9(a c (f3), (3) = for /3 > N. By observing that for every e > 0, 
P(i+eWi+e)j8 is stochastically dominated by P a ,p, we deduce that at c (/3(l+e)) > (l+e)a c (/3). 
Since by Theorem [TJ a c (/3) > for (3 G (N,N 2 ), this gives that a c (/3) is strictly increasing 
on (N,N 2 ). 

In order to prove continuity from the right, we note that for all 5 > 0, 9{a c {f3) + S, (3) > 
by definition. By the continuity of 9 (a, (3) we obtain that there exist e > 0, such that 
9{a c {0) + 5, /3 + e) > and therefore a c (/3 + e) < a c (/3) + 5. This together with a c (/3 + e) > 
a c ((3) completes the proof. □ 

Lemma 13 a c {0) is continuous from the left for (3 G (0,iV 2 ). 

To prove this, we use pQ. In that paper it is shown that for long-range percolation on Z d , 

inf{a : 9 (a, (3) > 0} = sup{a : E a<(3 (\C(0)\) < oo}. (16) 

Inspection of the proof of this result yields that this proof also works on the hierarchical 
lattice. Now use the following lemma. 

Lemma 14 Let a > and (3 > N . IfE a> p(\C(0)\) < oo, then there exist e > such that 
E Qi/3(1 _ £) (|C(0)|) < oo. 

Proof of Lemma [731 (given Lemma\L$. For every e > 0, we have that 

F M/3 ), /3 (i- £ )(|C(0)| = oo)>0, 

by the strict increase of a c ((3). This implies that R ac m\ph_ s \(\C(0)\) = oo, for every e > 
and therefore E Q , c ( /3 ) i/3 (|C(0)|) = oo. 

Furthermore, by fflBT) we know that for every 5 > 0, E Q!c ( / g)_5 ii g(|C(0)|) < oo. Therefore, 
there exist e > such that Mt ac ^_^ t pn^(\C(0)\) < oo, which implies that for all 5 > 0, 
there exist e > such that a c ((3(l — e)) > a c (f3)—5. This together with a c (/3(l —e)) < a c ((3) 
gives continuity from the left of a c {(3). □ 



Proof of Lemma 14_. Assign independent uniform(0, 1) random variables to all pairs of vertices 
in Qn. The random variable assigned to the pair (x,y) is denoted by U(x,y). We say that 
x and y share an edge for the parameters a and /3 if U(x,y) < 1 — exp(— a/3~ d ^ x ' y '). This 
construction provides a coupling for long-range percolation models with different values of 
a and (3. Define C(x;a, 0) as the cluster of vertices that can be reached by paths if the 
parameters are a and (3. 

Assume that a := E Q , )i g(|C(0)|) < oo and take e > small enough (we will see later 
exactly how small). Define *4o(0) := C(0;a, /?)■ In an inductive fashion, let *4.- +1 (0) be the 
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set of vertices not in U* =0 A, (0) that can be reached from Ai(0) by crossing an edge present 
for the parameters a and /3(1 — e). Note that ^(0) C Ai(0) and that by construction, 
C(0;ct,/3(1 - e)) = U^A(O). The next step in the proof is to bound E(|A(0)|). By 
definition E(|^4 (0)|) = a. Since the graph is transitive, for every x G Ai(0) the expected 
size of the set 

{yeSl N \ u; =0 A,(0); U(x, y)<l- exp(-a((l - 
is bounded above by 

b := p (f / (0, V) < 1 - exp(-a((l - e)f3)- d(0 ^) U(0, y) > 1 - exp[-ap- d(0 ' y) ] 

P(l - exp(-a/3- d ^) < U(0, y) < 1 - exp(-a((l - e)/3)- d(0 ' y) ) 
\ P(l - exp{-a/3- d (°^) < U{0,y)) 



VfiV - ivy»-i ex P(~ a ft ')-exp(-Qi^ '(1 - e) 
^ l J exp(-a/3- J ) 



i=i 

oo 



J2(N - 1)^(1 - exp(-«/3-((l - e)"* - 1))) 



i=l 
oo 



< J^(JV - l)^-^-^! - e)-* - 1] 



i=l 



(JV-l)a 7 ^./^ (iV-l)a 1 



l-e)0-N " /3 — N 

ae(N-l)N 



{/3{l-e)-N)W-N)' 

which converges to 0, if e \ 0. Therefore, we can choose e > so small that b < a -1 . Note 
that E(|^ +1 (0)|) < 6E(|A(0)|), and because of the transitivity of the graph, E(|A+i(0)|) < 
aE(|^ +1 (0)|). So, we have 

E(|A(0)|)<(a&yE(|A)(0)|). 
Since C(0; a,/3(l— s)) = U^ 1 >^(0) and ab < 1, we have 



C(0; a, p(l - e)) < V a(afeV = — ^— < oo. 

^— ' 1 — ao 



i=0 

□ 



Proof of Theorem^ The only things left to prove is that a c {(5) /* oo for (3 /* N 2 , but this 
follows immediately from Lemma HU equality ffl6|) and observing that a c (N 2 ) = oo. □ 
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7 Possible generalizations 



Possible generalizations of the model considered in this paper include: 

1. Randomness in the hierarchical lattice. The hierarchical lattice Qn is generated by 
a iV-regular tree. An interesting question is how randomness in the underlying tree 
(or induced random metric) affects the percolation process on the resulting lattice. 
One possibility is that the metric generating tree, is a Galton- Watson tree. Analysis of 
long-range percolation on such random hierarchical structures is not a trivial extension 
of the analysis in this paper, since in renormalisation schemes, one has to take care of 
all kinds of dependencies of the sizes of balls of given diameter. 



2. More general connection function p(k) . In this paper we focused on p^ = 1— exp 




What are necessary and sufficient conditions on g(k) so that whenp^ = 1— exp (-ag(k)) 
we have < a c < oo? 

3. Random cluster models. We only consider independent percolation on the hierarchical 
lattice. We did not try to incorporate Random cluster (or Fortuin-Kasteleyn) model 
[H] yet. Some work has already been done for the Ising model on the hierarchical 
lattice [T5] . 
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